Topological aspects of an exactly solvable spin chain 
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We analyse a spin-1/2 chain with two-spin interactions which shown to exactly solvable by Lieb, 
Schultz and Mattigi. We show that the model can be viewed as a generalised Kitaev model that 
is analytically solvable for all defect sectors. We present an alternate proof that the defect free 
sector is the ground state, which is valid for a larger parameter range. We show that the defect 
sectors have degenerate ground states corresponding to unpaired Majorana fermion modes and that 
the degeneracy is topologically protected against disorder in the spin-spin couplings. The unpaired 
Majorana fermions can be manipulated by tuning the model parameters and can hence be used for 
topological quantum computation. 
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I. INTRODUCTION 

There is much interest in the physical realisation of 
model systems for topological quantum computation 
(TQC)^. TQC is being actively explored as a physi- 
cal way to achieve fault tolerent quantum computation. 
In this scheme the qubits are realised using non-abelian 
anyons. The braiding operations on them implement 
quantum gates robustly. One of the simplest class of 
non-abelian anyons are realised in systems with unpaired 
Majorana fermions (UMF)'^. Kitaev presented a remark- 
able solvable spin-1/2 model on a honeycomb lattice^ 
which realises non-abelian anyons made up of UMF. The 
model can be written in terms of Majorana fermions in 
the background of Z2 gauge fields and the problem re- 
duces to solving a theory of non-interacting Majorana 
fermions in the background of static Z2 gauge field con- 
figurations. Kitaev showed that the ground state is the 
flux free configuration. The model has a phase which is 
characterised by a topological invariant, the Chern num- 
ber of the fermions, being equal to ±1. In this phase, 
there are UMF trapped to each vortex (a Z2 flux). Thus, 
if vortices can be created and manipulated, it is possible 
to braid the UMF. 

While anyons (abelian and non-abelian) are intrinsi- 
cally two dimensional objects, Alicea et. al showed that 
the braiding operation of UMF can be realised in one- 
dimensional wire networks^. They proposed a realisa- 
tion of such networks in semiconductor wires which can 
be engineered into a state supporting UMF at the edges. 

Kitaev's honeycomb model can easily be generalised 
to any lattice with coordination number three, if all the 
bonds can be coloured using three colours^rii. There are 
theoretical proposals to realise such systems in cold atom 
system&iS and Josephson junction quantum circuitsJ^. 
We had proposed and analysed such a generalised one 
dimensional model, which we called the Tetrahedral 
modeli^, where we showed that there are UMF trapped 
to defects and that the defects could be created and ma- 
nipulated by tuning the hamiltonian parameters. How- 
ever, this required the engineering of 3-spin interactions 



which is not easy experimentally. 

In this paper we study a simpler quantum spin-1/2 
chain which involves only 2-spin interactions and can 
thus, in principle, be realised in quantum circuits^'*. 
The model has been studied earlier by Leib, Shultz and 
Mattisi, who showed that the model was exactly solvable 
using the Jordan- Wigncr transformation. They showed 
that the model has infinite conserved 2-spin operators 
and proved that the ground state was in the (suitably de- 
fined) defect free sector. In each defect configuration, the 
model reduces to an analytically solvable fermion hop- 
ping problem and is characterised by degenerate ground 
states, the exact degeneracy depending on the defect con- 
figuration. 

We solve the model exactly using Kitaev's method**. 
We reproduce the previous results^ and extend the valid- 
ity of the regime of the proof of the ground state being the 
defect free sector. We then show that there are localised 
UMF zero energy modes in the sectors with defects. The 
ground state degeneracy corresponds to these modes be- 
ing occupied or unoccupied. The topological nature of 
these zero-energy modes makes them robust against dis- 
order in the strengths of couplings in the Hamiltonian. 
By tuning the couplings of the conserved 2-spin opera- 
tors in the Hamiltonian any defect sector can be made 
the ground state. Thus it is possible to manipulate the 
UMF and move them along the chain. 

The rest of this paper is organised as follows. In sec- 
tion ini we present the model, its conserved quantities 
and the Jordan- Wigner transformation which enables us 
to rewrite the theory in terms of Majorana fermions hop- 
ping in the background of static Z2 gauge field configu- 
rations. The diagonalisation of the Majorana fermion 
problem is described in section IIIII Section IIVI gives an 
exact analytic solution of the model. In section |Vl we 
present a proof that the groundstate lies in the defect 
free sector. Our proof has a larger regime of validity 
than the one given previously*. It also yields an expres- 
sion for fermionic gap. Section I VII contains a detailed 
analysis of the zero modes of the Majorana fermions. We 
show that the zero mode and, therefore, degeneracy in 
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the multiparticle spectrum is topological in origin and is 
robust under the variation of the couplings in the hamil- 
tonian. We summarise our results and discuss them in 
the concluding section IVlIl 



II. THE MODEL 



where Ui = -ir]i^ir]i^2 and (*?i,iCi,2) Ilf = are 
conserved quantities. 

The Majorana operators and rji follow the anti- 
commutation relations, 

{(i,^]} = 2(5,j;{77i,?7j} = 25ij;{^,,r]j} = 0. (6) 

Thus, the spin Hamiltonian gets converted into the 
fermionic Hamiltonian with periodic boundary condition 
when = +1 and with anti-periodic boundary condi- 
tions when = —1. 



FIG. 1. The XYZ-Ising chain. There are two sites per unit 
cell. The x, y and z bonds are as indicated. 



We consider the Hamiltonian, 

N 

H = J2 ('^^<1<2 + Mll^h + •^^<1<2 + <2<+l,l) 
i=l 

(1) 



III. DIAGONALISATION 

The Hamiltonian can be diagonalised in the standard 
way. We write the eigenstates as direct products of states 
\Q) in the 77 fermion sector and states \M.) in the ^ 
fermion sector. We will refer to states belonging to r] 
fermion sector as the gauge sector and states belonging 
to ^ fermion sector as matter sector. We choose the states 
in the gauge sector to be the simultaneous eigenstates of 
the Z2 flux operators, i.e \Q) = \{ui\), where 



with periodic boundary conditions. 



(2) 



Ui\{Ui}) = Ui\{Ui}). 



(7) 



In each unit cell, there is a conserved Z2 invariant. We then have 



Wi = erf ^cf 2- Apart from these N local conserved quan- 
tities, there are also three global conserved quantities cor- 
responding to global TT rotations about each of the three 
axes which are symmetries of the model. We denote these 
by, S° = rii^i f f,if f,2 These are not independent, we 
have = n»=i and S^S?' = 

Thus, we see that the Jz terms in the Hamiltonian ([T]) 
are invariants and, hence, do not affect the eigenstates. 
In this work we concentrate on = Jy = ^ ■ As we will 
show, at this point the model is analytically solvable. 
Also, note that this point the model has a global U{1) 
symmetry corresponding to spin rotations about the z 
axis. Leib et al had studied the model at the Jz ~ Jx ^ 
Jy point and had called it the Heisenberg-Ising model. 
We follow their nomenclature and refer to the model in 
equation ([!]) with — Jy as the XF-Ising model. 

We express the Hamiltonian in terms of Majorana 
fermions using the Jordan- Wigner transformation. 



where fi, = Ut<j '^li'^l2- 

Then, the Hamiltonian reduces to, 

X! ( ( 9^ ) *^^1^'.2 + JzU, 



(3) 
(4) 



(5) 



H{u,\ \M)\{ui)) =H{u,\ \M)\{ui)). (8) 
The single particle eigenvalue equation is given by. 



T-<Pi^i,2 H 7: <Pi,2 = e<pis 



-iJ(l - Mj) 



4 + 2'^»+i'i ^ <^'Pi,2^ 



with boundary condition 



'14 



(9) 



(10) 



We define the — 1 V i sector as the then defect 

free sector. This zero defect sector consists of two sectors 
corresponding to periodic or anti-periodic boundary con- 
ditions for Tlx = +^ and E^, = — 1 respectively. We will 
prove that this is the ground state sector of the model. If 
n of the Ui's are equal to -1-1 then we call it the n-defect 
sector. An n defect sector reduces to solving the hopping 
problem on n decoupled open chains. If we can solve 
one defect sector (one open chain), then we can solve all 
the sectors of the model because all of them reduce to 
decoupled open chain Hamiltonians. 



IV. EXACT SOLUTION IN ALL SECTORS 

In this section we put Jz — ^ since the eigenstates are 
independent of Jz- We first consider the zero-defect sec- 
tor of a chain with N unit cells, in the thermodynamic 
limit N ^ 00. The Hamitonian is diagonalised in terms 
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of the Fourier components of the Majorana fermion op- 
erators, 



fee 



(11) 



where Ri — ai, a being the lattice spacing, we have su- 
pressed the sub-lattice index and represent two compo- 
nent column vectors. The hermiticity of implies that 
= ^jj, their components satisfy the canonical anti- 
commutation relations, 



(12) 



Fermionic modes are thus defined on half the Brillouin 
zone. The momentum space Hamiltonian is. 



H 



1 



r/a 



dk ^Ih^k 



where/i(fc) = [ j _ ^ ^ 



ika 



(13) 
(14) 



The Hamiltonian in the diagonalised form is expressed as 
2 r^/^ 



H = - 



dk e{k) {ixtVxt - iXkVXk) (15) 



where, 



tanafc = 



-( ' 

— sin(fca) 



J — cos{ka) 
and e(fc) = - 2Jcos(fca) 4- 1. 



(16) 

(17) 
(18) 



Thus, the sytem is gapped for J ^ 1. The values of k are 
determined by applying the boundary condition (jlOp . 



ka — 



{2m — p)7r 
N ' 



(19) 



where m = 1, 2, 3 ... and p = 0, 1 for (PBC) and anti- 
(PBC) respectively. 

We now consider the one defect sector. The fermionic 
system is then an open chain with N unit cells satisfying 
the boundary conditions. 



and (l)N.2 — 0. 



(20) 



We use the standing waves to diagonalise this open chain 
single particle Hamiltonian, 



ikRi 



'2N 



')-k.ae 



-ikRi 



The boundary conditions ([201) imply that 



gjofc _ gi(-/V+l)fea 



(21) 
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FIG. 2. The red line shows the plot between k and n without 
the correction a^. The green and blue line shows the plot 
between k and n for J < 1 and for J > 1 respectively. 



Thus, k are determined by the solution of, 

riTT Oik 



ka 



N + 1 N + 1 



(22) 



We have solved equation (|22p numerically for = 10 
and plotted k versus n for J > 1 and J < 1 as shown in 
figure H 

We see the number of allowed values of fc is TV for J > 1 
but is iV — 1 when J < 1. As we will show later, the 
missing mode is a zero energy eigenfunction of the single 
particle Hamiltonian with the wavefunction peaked at the 
edges of the open chain. The relation J — e**^" = efce'"'= 
has been shown in the Argand plane in the figure ([3]) and 
figure (|H). The appearance of the zero mode depends 
on the topology of parameters of Hamiltonian as shown 
in figure ([3]) and (Hj). Zero mode appears when the unit 
circle shown there encloses the origin. This argument 
can be generalised to show that if the path circles the 
origin v times, then there will be v only N — v standing 
waves. V can be identified with a topological invariant of 
the closed chain in terms of the Berry potential defined 
on the half-Brillioun zone of the closed chain, 



A{k) = i {4>\dk(l)k - h.c 



(23) 



where (f>k are the two-component single particle wave- 
functions. The topological invariant is the winding num- 
ber of the relative phase of the two components of the 
wavefunction along the half Brillioun zone. 
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TT JO 



TT /a 



dk A{k) 



(24) 



The geometric interpretation of v is as follows. Con- 
sider a general Hamiltonian of the form 

h{k) = e{k)n{k) 
and n{k) = cos 0fcr^ -|- sin (cos fifcT^ -f sin f2fcr*')(25) 
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FIG. 3. J - e''"^ = efce^"'= in Argand plane. For J > 1, 
a(0) = and a[-K) = 0. 

Imaginary Axis 




Real Axis 



FIG. 4. J - e**^" = eke'- 
a{Q) = and a(Ti) — n. 



in Argand plane. For J < 1, 



where t° are the Pauli matrices. The Hamiltonian in 
equation p4|) is the case of dk — tt/2 and ilk — ctk- For 
fully connected near neighbour chains, the off-diagonal 
terms in the single particle Hamiltonian in equation do 
not vanish implying that 0^ 7^ 0,7r for any k. Further, 
if Efc 7^ the Hamitonian in equation (|25p represents a 
mapping of the half-Brillioun zone to the sphere with 
two points removed (the north and south poles in our 
convention), which is topologically equivalent to an an- 
nulus. Thus, v which is the winding number of the map is 
a topological invariant for the class of gapped chains for 
which no nearest neighbour hopping amplitude vanishes. 
Thus, we expect the zero energy edge state to be ro- 



bust and independent of disorder in the couplings, pro- 
vided the system remains gapped and fully connected by 
nearest neighbour couplings. We will demonstrate this 
explicitly in a later section. 

The solution of the open chain presented above can be 
used to to solve n defect sector, for any n as these sectors 
consist of n decoupled open chains. 



V. GROUND STATE AND GAP 

In this section, we will prove that zero defect sector is 
the ground state sector for = 0. In order to prove the 
statement we first prove that the ground state energy of 
zero defect sector is less than that of one defect sector for 
large N . We then prove that ground state energy of one 
defect sector is less than or equal to that of the n-defect 
sector. We then turn on Jz and compute regime of the 
stability of the system to defect formation. We finally 
compare our results with those of Lieb et. alX and point 
out precisely how our proof extends the regime of validity 
of their result. 

First step: We outline the proof that the ground state 
energy of zero defect sector GSE^ds is less that of the one 
defect sector, GSEods- The details are given in Appendix 
Rl The ground state energies of the two sectors are. 



GSEz 



GSEo 



N 

fc=i 

N 

fc=i 



2fc7r 

fcTT 

ivTT 



a{kTi/N + 1] 
N + 1 



(26) 



(27) 



In the thermodynamic limit of iV — >■ (X), we convert the 
summations into integrals using the Euler-Maclaurin for- 
mula to get. 



GSEzds = / e(fc)* 



GSEods — — 



SiJ) 



dk- 



e{k)dk + ,7(1 - S{J)) 

cos{k) — J 
^J2 - 2Jcos(fc) + 1 



(28) 



(29) 



We then show in Appendix that 6{J) < 1 for all finite J, 
thus proving that the ground state energy of zero defect 
sector is less than that of the one defect sector. 

Second step: Now, we show that the ground state 
energy of one defect sector is less than that of two defect 
sector and so on. To show this, we split up Hamiltonian 
into three parts. 



L1+L2 



H 



Li 



Hr 



Hr 



(30) 



where H is the Hamiltonian for one defect sector, H^^ 
and ^-re Hamiltonians of lengths Li,L2- H12 is the 
link bond between Hj^^ and H^^. 
From variational principle, 

eGiLi+L2) <{ib\HL,+HL,+Hi2W (31) 
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where ec [Li + L2) is the ground state energy of Hamil- 
tonian for one defect sector. 

Using the trial wave function ji/;) = iV'-Li) 1^1,2)1 g^t 

ec {L1+L2) <eG (ii) + eG (^2) 

+ {^lA{^l,\Hi2Hl,Ml,) (32) 

Let us consider the hnk Hamihonian H12 = 
where ^„ belongs to the Hamiltonian H^^ and ^„+i be- 
longs to the Hamiltonian Hj^^ . Then, 

{iPlA{'^L2\Hi2\^Li)\iI^L2) = i(V'Li|fn|V'Li)(V'L2lCn+l|V'L2)- 

Since the expectation value of single ^„ operator in the 
ground state is zero. Therefore, 

ea {Li + L2) < ea (Li) + EG (12) (33) 

Thus, we have proved that the ground state energy of 
one defect sector is less than the ground state energy of 
two defect sector. Similarly, we can prove for two defect 
sector and so on. Therefore, zero defect sector is the 
ground state sector at Jz = 0. 

We now turn on Jz and examine the region of stability 
of the zero defect sector under defect production. The 
difference of the ground state energies of the zero and 
one defect sector for non-zero Jz is 

GSEzds - GSEods = J(l - 5{J)) + 2Jz (34) 

Thus, the system will become unstable towards defect 
production when 

Jz<~^{l-5{J))- (35) 

This extends the result of Lieb et. al. since their proof 
is not valid for negative values of Jz- 



A. Nature of low energy excitations 

For Jz > --1(1 - 5{ J)), the zero defect is the ground 
state sector. The first excited state can be either the first 
excitation in the zero defect sector or the ground state 
of the one defect sector. The excitation gap in the zero 
defect sector follows from the single particle spectrum 
and is given by 

A = |(l-J)|. (36) 

Therefore, model is gapless for J — 1. Numerically, we 
have shown that for = first excited state of the 
model in the region < J < 0.75 and J > 1.17 is the 
ground state of the one defect sector . The excited state 
of the zero defect sector is the first excited state in the 
region 0.74 < J < 1.17. These results have been shown 
in figure ([5]) . We also see that the model has a finite gap 
for all non-zero values of J 7^ 1,0. 
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J 

FIG. 5. Low energy excitation in XY-Ising Model. The 
ground state energy of one defect sector is shown in green and 
first excitation of zero defect sector is shown in red. While 
plotting all the energies, ground state of zero defect sector has 
been taken as reference. 

VI. ZERO MODES 

We now analyse the zero mode of the Majorana 
fermions. When e = 0, equations ([9]) decouple and be- 
come 

^ + (Pi+1,1 = 

and - 0,_i^2 + ^^ii^lM^^ 2 ^ 0. (37) 
These recursion relations can formally be solved. 




There are two degenerate and independent solutions for 
every set of values of the parameters and every fiux con- 
figuration. In zero defect sector, imposing boundary con- 
dition (jji^i = (jjN+iA, we find that the zero modes exist 
for zero defect sector only at the gapless point J = 1. In 
one defect sector, which is an open chain with boundary 
condition 0o.i = (t>N.2 = 0, and 0^.2 can be expressed 
as, 

= ijy4>o,i 

and (j)i^2 = 4'N.2- (39) 

We can see that if J > 1, then (/)i_2 becomes zero at 
the boundary for infinite open chain. Now, if we take 
4'i,i = 0^*7 then the boundary condition at the other end 
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4>o,i = is also satisfied for ui = —1. Therefore, the only 
one solution for zero mode exist in an infinite chain for 
J > 1. The same is true for the case J < 1. In both 
cases, we get only one non-zero solution for zero mode. 

Thus, the zero mode occurs when the toplogical invari- 
ant is non-zero. It can also be seen from the recursive 
solution presented above that the zero modes exist even 
if the coupling, J varies with i. Thus the zero mode topo- 
logically protected against random, spatial variations in 
the exchange coupling. 



VII. CONCLUSION 

To summarize, we have analysed an exactly solvable 
spin- 1/2 chain which was earlier studied by Lieb, Shultz 
and Mattisi. We have shown that it is a generalised 
Kitaev chain. Using the methods of Kitaev^ we can re- 
produce the earlier results. The formalism also helps us 
extend the earlier proof to negative vlaues of Jz and com- 
pute the value of Jz when the defect free state becomes 
unstable to defect production. In addition we have ex- 
plicitly displayed the topological nature of the zero en- 
ergy modes in the presence of defects. We have identified 
the topological invariant that determines the existence 
of zero energy edge modes in the system. Thus the de- 
generacy is protected against random variations in the 
exchange coupling. 

In conclusion, our results show that zero energy Majo- 
rana modes can be created and manipulated in this model 
exactly in the same way shown in previous worki^. The 
system has a gap for all non-zero values of J 1. Thus 
the UMF can used as a topological qubit. The impor- 
tance of this model is that it has only two-spin couplings 
which can, in principle, be realised experimentall}»i^. 
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Appendix A: Details of proof 

The ground state energy GSEzds of zero defect sector 
is expressed as 

N 



GSEzds — — ^ e 



k=l 



2kn 



(Al) 



where e{k) = \J J"^ — 2Jcos{ka) + 1. 

To evaluate the summation series we convert it into 
integral using Euler-Maclaurin formula, 

Y^fix)^ / f{x)dx + B,[f{0)~f{n)] 

x=l 



where Bi = -1/2,^2 = 1/6,^3 = 0,^4 = -1/30, B5 = 
0,Be = 1/42... are the Bernoulli numbers. 

The expression for ground state energy for zero defect 
sector equation (IA1[) in thermodynamic limit N 00 
becomes. 



GSEzds / e{v)dv. 

1" Jo 



(A3) 



where we have used the identity e(0) — e(27r) — J + Jz. 
The ground state energy GSEods of one defect sector 

is. 



N 



GSEods — — ^ e 



fcvr a{k7r/N + l) 



N + 1 



N + 1 



(A4) 



To evaluate the summation series, we convert it into 
integral using Euler-Maclaurin formula and neglect last 
terms containing (N-|-l) in the denominator in limit N — 

00, 



GSEods 



N 







kn O'iwTi) 
N + 1 N+1 



dk 



1 

2^ \N + 1 



«(0) 



Ntt aj^,) 
N+1 N + 1 



(A5) 



In the thermodynamic limit N ^ 00, above equation 
[5]) becomes, 



( kn 



a{kiT/N + 1) 
iV+ 1 



dk 



(A6) 



Let us substitute y = ^"^^ , and above expression for 
ground state energy for one defect sector becomes 



GSEods^^^l ^[y + J^,)dy 

+ i[6(0)-eW]. (A7) 

Using Taylor expansion and neglecting terms containing 
1 

JV-l-l ' 



N+1 fTf+T 

GSEods = / <^{y)dy 

1" Jo 



1 f ™+i 



{a{y)e'{y))dy. (A8) 



Using integration by parts in the second term and then 
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taking the thermodynamic hmit, we get 

N r 



GSEods 



1 r , Jcos(2/)-j2 



^J2 - 2Jcos(y) + 1 



1 



+ - [a(7r)e(7r) +a(0)e(0)] 

TT 

From figures dSllll), 

e(0) = J + 1 

e(7r) = J — 1 when J > I 
e(7r) = 1 — J when J < 1 
a(0) = 

a(7r) = when J > 1 
a(7r) = TT when J < 1. 



(A9) 



Substituting all these values, we show that for both cases 
J > 1 and J < 1 ground state energy of one defect sector 
is expressed as, 



GSEods = GSE.ds + J(l - 5{J)) 
where 5[J) = - r dy^^M=^=.{An) 



7o - 2Jcos(y) + 1 



Since in the expression for S{J) the integrand is always 
less than 1, we can write 



5{J) < - [ dy=l. (A12) 
7^ Jo 



(AlO) 



Therefore, from equation (jAllI) . we can say that the 
ground state energy of zero defect sector is less than 
ground state energy of one defect sector. 
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